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ABSTRACT 

Estimates for the nonlinear a effect in helical turbulence with an applied magnetic field are 
presented using two different approaches: the imposed-field method where the electromotive 
force owing to the applied field is used, and the test-field method where separate evolution 
equations are solved for a set of different test fields. Both approaches agree for stronger fields, 
but there are apparent discrepancies for weaker fields that can be explained by the influence 
of dynamo-generated magnetic fields on the scale of the domain that are referred to as meso- 
scale magnetic fields. Examples are discussed where these meso-scale fields can lead to both 
drastically overestimated and underestimated values of a compared with the kinematic case. It 
is demonstrated that the kinematic value can be recovered by resetting the fluctuating magnetic 
field to zero in regular time intervals. It is concluded that this is the preferred technique both 
for the imposed-field and the test-field methods. 

Key words: magnetic fields — MHD — hydrodynamics - turbulence 



1 INTRODUCTION 



The a effect is commonly used to describe the evolutio n of the 
large-scale magnet i c field in hydromagn etic dynamos jMoffatd 
1 19781 : |Parkei)ll979l : iKrause & Radlei|[T98Q) . However, the a ef- 
fect is not the only known mechanism for explaining the gener- 
ation of large-scale magnetic fields. Two more effects have been 
discussed in cases when there is sh ear in the system: the incoher- 
ent alpha-shear dynam o t Vishniac & Brandenburg|IT997l : ISokolovl 
1 1 9971: ISilant' evi l200d: IProctor. ,20071) a nd the shear-current effect 



jRogachevskii & Kleeoriij|2003ll2004t) . In order to provide some 
understanding of the magnetic field generation in astrophysical 
bodies such as the Sun or the Galaxy, or at least in numerical sim- 
ulations of these systems, it is of interest to be able to identify the 
underlying mechanism. 

Astrophysical dynamos are usually confined to finite domains 
harboring turbulent fluid motion. Both the Sun and the Galaxy are 
gravitationally stratified and rotating, which makes the turbulence 
non-mirror symmetric, thus leading to an a effect. In addition, the 
rotation is nonuniform, which leads to a strong amplification of the 
magnetic field in the toroidal direction, as well as other effects such 
as those mentioned above. Instead of simulating such systems with 
all their ingredients, it is useful to simplify the setup by restricting 
oneself to Cartesian domains that can be thought to represent a part 
of the full domain. At low magnetic Reynolds numbers, i.e. when 
the effects of induction are comparable to those of magnetic diffu- 
sion, the a effect can clearly be identified in simulations o f convec- 
tion in Cartesian domains: see Brandenb urg et al. Il ll990l) . Here, a 
has been determined by applying a uniform magnetic field across 
the simulation domain and measuring the resulting electromotive 



force. This method is referred to as the imposed-field method. How- 
ever, in subsequent years simulations at larger magnetic Reynolds 
numbers have revealed problems in that the resulting a becomes 
smaller and strongly fluctuating in time. This was first found in sim- 
ulations wh ere the turbulen ce is caused by an externally imposed 
body force jCattaneo & Hu ghes 1996), but it was later also found 
for convection jCattaneo & H ughes 2006). This suggested that the 
mean-field approach may be seriously flawed (Cattaneo & HughesI 
I2OO9I) . 

Meanwhile, there have been a number of simulations of con- 
vection where large-scale magnetic fields are being generated. 
Such systems include not only simula tions in spherical shells 
te rown ing et alj2006l:lBrown et al .l2007h. but also in Cartesian do- 
mains (Kiipvla et al. 2008, 2009^: iHughes & Proctoill2009l) . How- 
ever, the absence of a significant a effect in some of these simula- 
tions led Hughes & Proctor ( 2009) to the conclusion that such mag- 
netic fields can only be explained by other mechanisms such as the 
incoherent alpha-shear dynamo or the shear-current effect. Such 
an explanation seems to be in conflict with earlier c laims of a finite 
Q effe ct as d etermined by the test-field method of ( Schrinn er et al.l 
20o3, 2007), and in particular with recent results for convection 



(Kiip vla et alj2009bh . The purpose of the present paper is therefore 
to discuss possible reasons for conflicting results that are based 
on different methods. The idea is to compare measurements of 
the a effect using both the imposed-field method and the test- 
field method. We consider here the case of helically forced tur- 
bulence in a triply-periodic domain. This case is believed to be 
well understood. We expect a to be catastrophically quenched, 
i.e. a is suppressed for field strengths exceeding the Zeldovic^ 
( Il957i) value of R^^^^Bcq, where Bcq is the equipartition field 
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strength where kinetic and magnetic energy densities are compara- 
ble JTiejmportance of the Zeldovich field strength was emphasized 
by iGruzinov & Diamond! ( Il994h in connection with catastrophic 
quenching resulting from magnetic helicity conservation. 

In this paper we focus on the case of moderate values of _Rni 
of around 30. This is small by comparison with astrophysical ap- 
plications, but it is large compared wi th the critical value for dy- 
namo action in fully helical turbulence dBrandenburdllOOlh . which 
occurs for JJm ^ 1 in our definition of _Rm based on the wavenum- 
ber of the scale of the energy-carrying eddies, i.e. the forcing 
wavenumber. In addition, we only consider cases with a magnetic 
Prandtl number of unity. However, this should not worry us too 
much, because we know that the large-scale dynamo wo rks inde- 
pendently of th e value of the magnetic Prandtl number jMininnil 
l2007l : lBrandenburg 2009) . 



2 HELICAL TURBULENCE AND a EFFECT 
2.1 Forced turbulence simulations 

Throughout this paper we consider hydromagnetic turbulence in the 
presence of a mean magnetic field Bq using triply-periodic bound- 
ary conditions. The total magnetic field is written as Bq + V x A, 
where A is the magnetic vector potential. We employ an isothermal 
equation of state where the pressure is proportional to the density, 
p = pel, with Cs being the isothermal sound speed. The govern- 
ing evolution equations for logarithmic density liip, velocity U, 
together with A, are given by 



Dlnp 
D[7 



Cs V In p, 



(1) 



(2) 



(3) 



— = J X (Bo + B)/p + f + Fvi. 
f)A 

— = Ux {Bo + B)+r,V^A, 
ot 

where Bo + B is the total magnetic field, but since Bo = const 
it does not enter in the mean current density, which is given by 
J — ^ X B / fiQ, where /xq is the vacuum permeability. Fur- 
thermore, D/Dt = + t/ • V is the advective derivative, 
f vise ~ p^^V ■ 2piyS is the viscous force, v is the kinematic 
viscosity, Sij — ^{Uij + Uj^i) — |(5ijV ■ U is the traceless rate 
of strain tensor, and / is a random forcing function consisting of 
plane transversal waves with random wavevectors k such that |fc 
lies in a band around a given forcing wavenumber kf. The vector 
k changes randomly fro m one timestep to the next. This method is 
described for example in 'Hansen et ail i2004h . The forcing ampli- 
tude is chosen so that the Mach number Ma = Urms/cs is about 
O.I. 

We consider a domain of size x Ly x L^- We use Lx = 
Ly = Lz = 27r/fci in all cases. Our model is characterized by the 
choice of magnetic Reynolds and Prandtl numbers, defined here via 



Ru 



(4) 



We start the simulations with zero initial magnetic field, so the field 
is entirely produced by the imposed field. The value of the magnetic 
field will be expressed in units of the equipartition value 



Bcq = {^lopu ) 



2\l/2 



(5) 



We consider values of Bq /-Bcq from 0.06 to 20 along with a mag- 
netic Reynolds number of about 26, adequate to support dynamo 
action. 



2.2 Q from the imposed-field method 

The present simulations allow us to determine directly the a ef- 
fect under the assumption that the relevant mean field is given by 
volume averages, denoted here by angular brackets. Given that the 
magnetic field is written as B = V x where A is also triply pe- 
riodic, we have (B) =0. We can determine the volume-averaged 
electromotive force. 



{£) = = {ux b), 



(6) 



where u = U (U) and 6 = B are the fluctuating components 
of velocity and magnetic field, and (B) = (V x A) = 0. 

For mean fields defined as volume averages, and because of 
periodic boundary conditions, we have (J) = 0. Under isotropic 
conditions there is therefore only the a effect connecting {£) with 
Bq via {£) = QimpBo, so 



Qimp = {£} ■ Bo/B'o. 



(7) 



In all cases reported below we assume Bo = (Bo, 0, 0). Note that 
V X (f ) = and therefore our time-constant imposed field is 
self-consistent. 



2.3 a from the test-field method 

A favored method of determining the full ag t ensor is by using the 
test-field method /Schrinner et alJ bOOl lioOTl) . where one solves, 
in addition to equations ([TJ-llSjl, a set of equations. In the special 
case of volume averages this set of equations simplifies to 

da" — 



^^Uxb'' + ux{Bo + B'')+uxb'' -u x ^r/V a", (8) 

where b" = ^ x a" with q = 1 or 2 denotes the response to each 
of the two test fields B". Throughout this paper, overbars denote 
planar averages. Later we consider arbitrary planar averages and 
denote their normals by superscripts, but here we restrict ourselves 
to xy averages. We use two different constant test fields, 



B^ = (S,0,0), B^ = (0,S,0), 



(9) 



where B = const is the magnitude of the test field, but its actual 
value is of no direct significance, because the B factor cancels in 
the calculation of a. 

However, given that the test-field equations are linear in b", 
this field can grow exponentially due to dynamo action. When 
becomes larger than about 20 times the value of B, the determina- 
tion of Q becomes increasingly i naccurate, so it is advisable to re- 
set b" to zero in regular intervals l lSuretalj200l) . We calculate the 
corresponding values of the electromotive force (f = {u x b") 
to determine the components 

Q,, = {£)l/B. (10) 

This corresponds to the special case fe = when consider- 
ing sinusoidal and cosinusoidal test func tions described elsewhere 
( [Brandenburg. Radler & Schrinneill2008h . 

Even though the test-field equations themselves are linear, the 
flow field is affected by the actual magnetic field (which is differ- 
ent from the test field), so the resulting a tensor is being affected 
("quenched") by the magnetic field. This was successfully demon- 
strated in lSrandenburg et alj ( l2008bl) , where Uij takes the form 



Oli 



a-i5ij + a2BiBj 



(11) 



Here B = B/|B| is the unit vector of the relevant mean mag- 
netic field. In the induction equation the a effect occurs only in the 
combination 
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Table 1. Overview of the different types of fields and tlieir meaning. 



field 



symbol magn induct, eqn test-field eqn 



imposed field 
meso-scale field 
test field 

test field response 



Bo 
B 



Bi 
B 



yes 
yes 



yes 



yes 
yes 



(12) 



and this is also what is determined by the imposed-fleld method, 
but it is different from the mean values of the components of the 
aij tensor. On the other hand, in the case of a passive vector field 
it is the mean components of Uij rather than the components of 
ajjBj that are of immediate importance jTilgner & Brandenburd 
I2OO8I) . 



2.4 a in the presence of meso-scale fields 

The relevant mean field may not just be the imposed field with 
wavenumber fc = 0, but it may well be a field with wavenumber 
fe = fci . Such a field would vanish under volume averaging, but it 
would still produce finite values of {BiBj). For the diagonal com- 




10.0 



Figure 1. Plot of the integrals and hiP)- 



Table 2. Summary of the expressions for tx (/3) and ty (fi) as well as e^, (0) 
and ey(O) for the X, Y, and Z branches. 



Branch 



e^(0) ey{0) 



ponents of {oiij) we can write 

{ctxx) = ai + €xa2, {ctyy) = ai + eya2, 



X p'^/ii + p'^) i/[2{l + /32) 
Y /i(/3) 



1/2 
1/2 
1/2 1/2 



where the factors 



and 



(13) 



(14) 



by integration. Thus, we write = /i (/3) where we have defined 

^'^^'-], /32 + 2/3sine + l'^^-\l-l/2/32 > 1, 

where 9 — ky has been introduced as dummy variable. Since By = 
in this case, we have ey — 0. 

Finally for the Z branch, where the wavevector of the Bel- 
trami field points in the z direction, we have B = (Bq + 
Bi cos kz, Bi sin kz, 0), we find ~ Ii{(3) and ey = I2{f3) with 



HP) = 



cos^6i 



d9 



/32 + 2/3 COS 61 + 1 27r 



Mf3) 



quantify the weight of the 02 term. For a purely uniform field 
pointing in the x direction we have ex — 1 and ey = 0, while 
for a Beltrami field of the form B — (cos kz, sin kz, 0) we have 
= = 1/2. 

In practice we will have a mixture between the imposed field 
(below sometimes referred to as large-scale field) and a dynamo- 
generated magnetic field with typical wavenumber fc = fci (below 
sometimes referred to as meso-scale magnetic field). The solution 
to the test-field equations, b'', can also develop meso-scale fields 
with wavevectors in the x 01 y directions, but not in the z direc- 
tion, because that component is removed by the term u x bi in 
equation ([8}. Table [T] highlights the difference between imposed, 
meso-scale, and test fields. We denote the ratio of the strengths of 
imposed and meso-scale fields as /3 = _Bo /Bi and distinguish three 
(and later four) different cases, depending on the direction of the 
wavevector of the Beltrami field. 

The first case is referred to as the X branch, because the 
wavevector of the Beltrami field points in the x direction. To cal- 
culate ex there is, in addition to the imposed field Bo, a Bel- 
trami field _Bi(0, cos fcx, sin fcx), which does not have a com- 
ponent in the x direction. Thus, Bx = Bq, and since B = 
(Bo, Bi cos kx, Bi sinkx), we have — Bq + Bi, so ex = 

Bl = Bl/{Bl + Bl), or ex = f3^/{l + 0^). Likewise, with 0.06 < i?m^^Bo/Beq < 20 for ^ 26 and Pm = 1. In all 



where /o(/3) = (l+/3^)/[2(l-/3^)] and 61 = fcz has been used as a 
dummy variable. A graphical representation of the integrals is given 
in Fig. [T] and a summary of the expressions for ex{P) and ey{l3) 
as well as exiQ) and ey(0) for the X, Y, and Z branches is given 
in Table 12] The singularity in /o(/3) could potentially affect ctyy. 
However, the results shown below show that, at least for stronger 
fields, a2 goes to zero near the singularity of Io{P) such that ayy 
remains finite. 



3 RESULTS 

We have performed simulations for values of Bq in the range 

1/2, 



By = Bl cos kx we find for the volume average or, in this case, 
the X average {Bl) = B'}/2, so ey = 1/[2(1 + /3^)]. 

The next case is referred to as the Y branch, because the 
wavevector of the Beltrami field points in the y direction. Thus, 
we have B = {Bo + B\ sin ky, 0, Bi cos ky), so B^ = Bq + 
2B0B1 sin ky + Bl.lhis. is no longer independent of position, so 
the volume average or, in this case, the y average has to be obtained 



cases we use kt/ki = 3, which is big enough to allow a meso-scale 
magnetic field of wavenumber k\ to develop within the domain; see 
Fig.|2] We did not initially anticipate the importance of the meso- 
scale fields. Different runs were found to exhibit rather different 
behavior which turned out to be related to their random positioning 
on different branches. We used the existing results from different 
branches as initial conditions for neighboring values of Bq. 
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X branch Y branch Z branch YZ branch 




Figure 2. Visualization of Bz on the periphery of the computational domain for the X branch and Bx for the Y, Z, and YZ branches. The coordinate directions 
are indicated on the first panel. 



In this paper, error bars are estimated from the averages ob- 
tained from any of three equally long subsections of the full time 
series. The error bars are comparable with the typical scatter of the 
data points, but they are not shown because they would make the 
figure harder to read. Note that the results in this section consider 
saturated fields. The opposite case will be considered in Sect.|4] 



3.1 Different branches 

The resulting values of a are shown in Fig.[3l For strong imposed 
magnetic fields, RmBo/ B^^ > 1, the resulting dependence of a 
on Bo obeys the standard catast rophic quenching formula for the 
case of a uniform magnetic field jVainshtein & Cattanedll99"3) . 



Qflt — 



l + Rn,ByB, 



(for B = So = const only), 



(15) 



where ao — —^Urms is the relevant kinematic reference value 
for fully helica l turbulence with negative helicity and i?ni > 1 
jSur et al]|2008l) . We treat Rm as an empirical fit parameter that 
is proportional to Rm and find that Rm/Rm ~ 0.4 gives a reason- 
ably good fit; see the dash-dotted line in Fig. |3] The existence of 
such an empirical factor might be related to fact that the relevant 
quantity could be the width of the magnetic inertial range, and that 
this is not precisely equal to Rm- For RmBo/ B^^ > 1, a similar 
result is also reproduced using the test-field method, although a^x 
is typically somewhat larger than ciimp- 

For weak imposed magnetic fields, RmBf)/ B^^ < 1, ap- 
parent discrepancies are found between the imposed-field method 
and the test-field method. In fact, in the graphical representation in 
Fig.[3]the results can be subdivided into four different branches that 
we refer to as branches X, Y, Z, and YZ. These names have to do 
with the orientation of a dynamo-generated magnetic field. These 
dynamo-generated magnetic fields take the form of Beltrami fields 
that vary in the x, y, and z directions for branches X, Y, and Z, 
while for branch YZ the field varies both in the y and z directions. 
Earlier work without imposed fields has shown that branch YZ can 
be accessed during intermediate times during the saturation of the 
dynamo, but it is not one of the ultimate stable branches X, Y, or Z. 

Branches Y and Z show the sudden onset of suppression of 
Qimp for weak magnetic fields. This has to do with the fact that for 
weak imposed magnetic fields a dynamo-generated field of Bel- 
trami type is being generated. Such fields quench the a effect, even 
though they do not contribute to the volume-averaged mean field. 
On branch YZ the a effect is only weakly suppressed, while on 
branch X the imposed-field aimp increases with decreasing values 
of Bo. 
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Figure 3. Volume-averaged values of aimp, ctxx, and ctyy. A tilde indi- 
cates that the values are normalized by OQ, i.e. CKimp — '-'^imp /^O (solid 
line), OLxx = {axx)/ao (dashed line), 

^yv — i'^yy) / ^0 (dotted line), 
and Qfit = Ofit / ao (thick gray line, but only shown in the second panel). 
The two open symbols in the top panel indicate that the values of olxxIolo 
are negative. 
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Figure 4. Root-mean-square values of the mean magnetic fields as func- 
tions of the imposed field for turbulence with Rm = 26 for the X, Y, Z, 
and YZ branches in the same order as in Fig. [3] Diamonds, triangles, and 
squares denote B^^\ B'''', and B^^\ respectively. 

The test-field method reveals that on branches X, Y, and 
YZ the ayy component is nearly independent of Bo, and always 
larger than the axx component. However, on branch Z and for 
JJmSo/^cq < 1 we find that a^x ~ ctyy and only weakly sup- 
pressed. 

A comment regarding the discontinuities in Fig. [3] near 
i?niSo/5cq = 1 is here in order. The systems considered here 
are in saturated states. To the left of the discontinuities the system 
has a saturated meso-scale dynamo, while to the right there is none. 
Intermediate states are simply not possible. Hence, the discontinu- 
ities are caused by the effects of the meso-scale magnetic fields on 
Wrms and thus on 7?m. 

3.2 Relation to ai and 0:2 

In the following we will try to interpret the results presented above 
in terms of equation illl and determine qi and 02 for the different 
branches. For small values of Bo, a magnetic field with k = ki and 



Figure 5. Dependence of ai and 02 on Bq for the X, Y, Z, and YZ 
branches in the same order as in Fig. [5] 

hence a finite planar average can develop. Compared with the large- 
scale field Bo, we refer to this dynam o-generated field as meso- 
scale magnetic field. As demonstrated in lBrandenburd ( l200ll) . three 
types of such mean fields are possible in the final saturated state. 
These fields correspond to Beltrami fields of the form 




where = cos(fci^ + 0) and = sin(fci^ + (j!)) denote cosine and 
sine functions as functions of ^ = a;, y, or z, with an arbitrary phase 
shift 00 The precise value of Bi emerges as a result of the simula- 
tion, but based on simulations in a periodic domain jBrandenburd 
I2OOII) we know that Bi/Bcq should be about (fef/fei)^/^ times the 
equipartition value. This is also confirmed by the present calcula- 
tions. 

^ Unlike the case considered by [Brandenburg et al.l j2008bl) , here the test 
field has k = 0, and there is no relative phase to be considered. 
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Let us now discuss separately the different branches. As can 
be seen from Fig.|4] the weak-field regime is characterized by the 
presence of meso-scale magnetic fields that vary either in the x 
direction (the X branch), the y direction (Y branch), the z direction 
(Z branch), or in both the y and z directions (YZ branch). 

In order to get some idea about the values ai and Q2 on the 
various branches, we consider two limiting cases. For strong im- 
posed fields, j3 —* 00, the results lie formally on the YZ branch 
branch, because such a field has only very little variation in the 
X direction. However, {BiBj} will be dominated only by the uni- 
form field in the x direction, so we have ex = 1 and = 0; see 
Sect. l2.3l This means that fiimp = a^x = cii + Q2 and ciyy = ai, 
so we can calculate 



Ql 



0L2 



(17) 



where a tilde indicates normalization by qq- For weak imposed 
fields, /3 ^ 0, we can calculate Si and 012 on the X branch by 
using using the relations 



— Ql + \0L2, 

— Ql + Q2- 



(18) 
(19) 
(20) 



However, on the X branch ctxx is ill-determined, as seen in Fig. [5] 
and discussed in Sect. 14.11 below. Therefore we use only equa- 
tions l |19t and OOl l to calculate 



Ql 



2a„ 



Q2 = 2ain 



(21) 



For the Y, Z, and YZ branches, on the other hand, these relations 
have to be substituted by 



Ql 



2Qa;j3 



OL2 = 2ainip — 2(5 1 



(22) 



The resulting values of ai and 02 are plotted in Fig. |5] for each 
of the four branches. On the Y branch one can, as a test, also use 
the independent relation ol\ = Q.yy. The resulting values are about 
50% larger than the values shown in Fig. |5] suggesting that there 
could be additional contributions in the simplified relation &yy = 
Ql. On the Z branch, of course, axx = oiyy, so here too we have to 
use the equations i22t . 

In all cases we find that a is quenched by ol\ and 0.2 hav- 
ing opposite signs and their moduli approaching each other This 
is particularly clear in the case of strong fields where qi and 
— Q2 become indistinguishable, while each of them is still increas- 
ing. We note that th e turbulence itself is not strongly affected 
( [Brandenburg & Subramanian 2005a) . On the Y and Z branches 
both Ql and 02 are of order unity, but on the X branch they can 
reach rather large values when the imposed field is weak. The be- 
havior on the YZ branch is somewhat unsystematic, suggesting that 
this branch is really just the result of a long-term tra nsient, as was 
alread y found in the absence of an imposed field jBrandenburd 
l200lh . However, we decided not to discard this branch, because 
it is likely that transient solutions on this branch may become even 
more long-lived as the magnetic Reynolds number is increased fur- 
ther. 



3.3 Enhancement of Qimp in the field-aligned case 

The suppression of q = qi -I- Q2 by the magnetic field is not sur- 
prising. What is unexpected, however, is the dramatic enhancement 
of both Ql and — Q2 for weak imposed fields and equipartition- 
strength meso-scale fields that vary in the x direction (the field- 
aligned case or X branch). In this case the interactions of the cur- 



rent density associated with the Beltrami field and the imposed field 
generate a force varying along x, perpendicular to the components 
of the meso-scale Beltrami field. This generates a meso-scale ve- 
locity that in turn damps the Beltrami field, resulting in the slower 
rise in S'^' as Bo/Sgq is decreased. Further, the cross-product 
of the meso-scale velocity field with the Beltrami field generates a 
large-scale electromotive force in the x direction. This is seen both 
in Qimp and in axx- A rough estimate of this electromotive force 
can be obtained by considering the fields 



Bo 




Bi 



/ ^ 

B\ cos kx 
\ sin kx , 



(23) 



so that /loJi = ~kBi, where subscripts 1 denote meso-scale 
fields. The meso-scale current density and the imposed field will 
generate a meso-scale Lorentz force which will drive a meso-scale 
velocity field Ui. We estimate JJi by balancing 



Ji X Bo/p + utV'Ui ^ 0, 



(24) 



where is the turbulent viscosity. We therefore expect that U 
will saturate for 

^* \ — cos kx / 



U 



(25) 



This velocity field will generate an £0 parallel to Bo in conjunction 
with Bi 

£0 = {Ul X Bl) = QmcsoBo, (26) 

with Qmeso = B\ / {p/iQi^tk). We then expect the total Qimp to be 



Qi 



Q + 



(27) 



Normalizing by qq = — Mrms/3 and assuming Ut « Urms/Sfcf we 
find for small imposed field and a meso-scale dynamo that varies 
along x: 



Qin 



Qo 



9 



ki 



Bj_ 

Bca 



(28) 



Given that fcf/fci — 3 and noting that Bi/B^q reaches values up 
to 1.2, we find that Qimp/Qo ~ 40, which is still somewhat below 
the actual value of 53, see the top panel of Fig. [3] The remaining 
discrepancy may be explicable by recalling that the actual value of 
ft may well be reduced due to the presence of an equipartition- 
strength magnetic field. 



3.4 Comment on wavenumber dependence 

In previous work on the test-field method we used test fields with 
wavenumbers different from zero. It turned out that in the kinematic 
regime, a is proportional to l/fl + a(fc/fcf )^], where a = .5, 1 
l lBrandenburg, Radler & SchrinneJlBosI ; iMitra et aljlioogi) . It was 
shown that the variation of a with k represents nonlocality in space. 
In order to get some idea about the dependence of axx and ayy on 
k in the present case we compare in Table [3] the results for k — ko 
with those for fc = 0. It turns out that both values decrease by 30% 
on the X branch, and increase by less than 10% on the Z branch. 

The k dependence for the Z branch is minor, although one 
would have expected a small decrease rather than an increase. Nev- 
ertheless, within error bars, this result is possibly still compatible 
with the dependence in the kinematic case. For the X branch the 
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Table 3. Examples of the dependence of a^x and ayy on the wavenumber 
k of the test field. Note that the field strength is different in both cases. 



Branch fe/fco 



Ctx 



0.72 ± 0.14 
0.61 ± 0.02 
0.34 ±0.02 
0.35 ± 0.01 



0.51 ±0.16 
0.37 ±0.01 
0.32 ±0.02 
0.35 ±0.02 



0.06 
0.06 
0.2 
0.2 



error bars for fc = are larger. This is because of the strong in- 
teraction between the imposed uniform field and a Beltrami field 
varying along the same direction, as discussed in Sect. 13.31 It is 
therefore not clear whether the k dependence is here significant 
and how to interpret it. 



4 RESETTING THE FLUCTUATIONS 
4.1 Effectiveness of resetting the fields 

The evolution equations used both in the imposed-field method 
and in the test-fie ld method allow for dynamo action. This led 
lOssendriiver et al.l (2002) and Kiipylaetal. (2006) to the tech- 
nique of resetting the resulting magnetic field in regular intervals. 
This method is now also routinely used in the test-field approach 
jSur et al]|2008l) . and we have also used it throughout this work. 
The lack of resetting the magnetic field may also be the main 
reason for the rather low v alues of a found in the recent work 
of jH ughes & Proctorl l200St) ; see the corresponding discussion in 
iKapvla et al.u2009bir 

In this section we employ the method of resetting B to obtain 
better estimates for a for weak imposed fields, and to compare this 
with results from the test-field method. The result is shown in Fig.|6] 
where we show the dependence of aimp on Bo and on the reset in- 
terval At. We note that, in units of the turnover time, the reset inter- 
val Atuimskf has a weak dependence both on Bo and At, because 
small values of Bo and At quench itrms only weakly. The resetting 
technique has eliminated the branching for weak fields. For weak 
fields we find that the value of Qimp is slightly below ao, but this is 
partly because for finite scale separation there is an additional fac- 
tor {l+kl/kfy^ ^ 0.9 (Brandenburg, Riidler & Schrinner 20Qi). 
The actual value of Ofimp is somewhat smaller still, which may be 
ascribed to other systematic effects. 

It turns out that over a wide range of reset intervals the result- 
ing values of Qimp are not dep endent in a systematic way on the 
reset interval (see also iMitra e t al. 2009), although it is clear that 
the error bars increase for larger values of At. The same is true for 
the values of a^x and Uyy obtained using the test-field method, ex- 
cept for the case of weak fields on the X branch where the values 
of Oxx are ill-determined; see Table |4l where we compare the val- 
ues of Uxx and ayy for two different reset times in the case where 
Oxx is found to change sign (i?m^Bo/Bcq ~ 0.2). The increas- 
ing fluctuations for longer reset intervals occur as the system exits 
the kinematic regime. It might therefore be possible to find indi- 
cators of when the kinematic regime has been exited and resetting 
becomes necessary. However, we have not pursued this further in 
this work. 

For even larger values of At there is enough time for the meso- 
scale magnetic field to develop. An example is shown in Fig. |7] 
where 18 intervals of length Atttmisfcf — 270 are shown. For half 



1.00 



0.10 



0.01 





10 



100 



Figure 6. Dependence of Qimp (solid lines) and ogt (dotted lines) on the 

imposed field strength with fixed reset time AtUrmski = 50, ...,70 (upper 

1/2 

panel) and the dependence of oimp on the reset time for i?,,; Bo / Boq = 



0.1 (lower panel). In all cases we have 



' 30. 




1000 2000 3000 



4000 



5000 



270 with 



Figure 7. Time series of CKinip for j\tU];mskf 
1/2 

Rm Bo/Bcq = 0.1. The reset intervals are indicated by dotted vertical 
lines. In all cases we have Rm ~ 30. 



of these intervals the wavevector of the Beltrami field begins to 
develop in the x direction, so aimp is heading toward the X branch. 
In the other half of these cases the magnetic field is weak and Qimp 
lies on one of the other branches. None of these cases reproduce the 
correct kinematic value of a, because we are not really considering 
a kinematic problem in this case. This underlines the importance of 
choosing reset intervals that are not too long. 
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Table 4. Comparison of the results for a^x and ayy for two different reset 

1/2 

times At for the examples of the X and Z branches with Bq / Bcq = 
0.2. The reset time is normalized by the inverse turnover time (ttrmsfcf )~^. 



Branch At Ur 



25 -0.08 ±0.13 0.54 ±0.02 

50 -0.98 ±0.09 0.70 ±0.04 

25 0.34 ± 0.02 0.32 ± 0.02 

50 0.32 ±0.01 0.33 ±0.03 



50 



100 




150 



200 



Our results support the hypothesis that the precise value of the 
reset time interval is not critical except for the field-aligned case 
where the diagonal components of the aij tensor are large and quite 
uncertain, as indicated also by the large error bars. The sign-change 
found for axx at low or intermediate field strengths might therefore 
not be real. 



4.2 Time averaging in the test-field method 

We have already demonstrated that the length of the reset interval 
is not critical for the value of a, but longer reset times tend to lead 
to larger errors. In the present section we demonstrate this for the 
test-field method using the idealized case where the turbulent flow 
velocity is replaced by simple stationary flow given by the equation 



U = kfipz + V X {ipz), 
with 

ip = 'fi{x, y) — uo cos kox cos koy, 



(29) 



(30) 



which is known as the Roberts flow. 

When the magnetic Reynolds number exceeds a certain crit- 
ical value of around 60, some kind of dynamo action of 6' 
conmiences. This typ e of dynamo is often referred to as small- 
scale dynamo action jSur et al.l l2008l : [Brandenburg et alj l2008bl : 
ICattaneo & Hughe j|2009l)~ but this name may not always be ac- 
curate. In the case of the Roberts flow there would be no such dy- 
namo action if the wavenumber of the test field is zero, A: = 0, as 
assumed here. However, for k — ko, for example, dynamo action 
for the test-field equation is possible. The test fields are therefore 
chosen to be 

' cos kz \ -^2 ( sin kz \ 



B 







B 



cos kz 



w 

B 



(32) 



■ \ 
sin kz I , 
^0 / - \ / 

see lSur et al Since now the mean fields are also functions 

of z, the term u x b i cannot be omitted in equation jS) . 

As stressed by [Brandenburg. Radler & Schrirmej ( |2008|) . in 
the expression for the electromotive force there is in general also 
a contribution £o that is independent of the mean field. Given that 
test fields i?* are independent of time, we have 



£''{z,t) = £l{z,t) + a{z)B''{z) - r|,{z)^loJ''{,z), 



(33) 



where overbars denote xy averages (not volume averages), so there 
is also a term rjtfioJ'', where rjt is the turbulent magnetic diffusiv- 
ity. We have assumed that a and rjt are independent of time, and in 
this case they are also independent of z. The £q{z, t) term can be 
eliminated by averaging over time, i.e. (f q) = 0, so 



1.0 



0.5 



s 0.0 
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Figure 8. Plot of the instantaneous a for Rni = 65 (upper panel) and 
Rm = 65 (lower panel). In both cases running means are overplotted and 
converge to nearly the same value of about —0.096 in the upper panel and 
—0.090 in the lower one. The envelope functions are well described by 
exponentials and are also overplotted. Note however the different scales on 
the ordinate of both panels. The dash-dotted line shows the zero level. 



(34) 



In Fig. [8] we show the evolution of a for the Roberts flow with 
i?m = 65 and 55. In the case with i?m = 65 there are exponen- 
tially growing oscillations corresponding to a wave traveling in the 
z direction. In general such fields can be a superposition of waves 
traveling in the positive and negative z directions. It is seen quite 
clearly that the running time average is stable and well defined. The 
results for i?ni = 65 and 55 are close together (a/ao = 0.096 and 
0.090, respectively), suggesting continuity across the point where 
dynamo action sets in. This supports the notion that averaging over 
time is a meaningful procedure. 



5 CONCLUSIONS 

The present simulations have shown that the imposed-field method 
leads to a number of interesting and unexpected results. For im- 

— 1/2 

posed fields exceeding the valu e flm J3eq one recovers the c atas- 
trophic quenching formula of I Vainshtein & Cattaneol jl992h ; see 
equation l |15t . We emphasize once more, however, that this for- 
mula is only valid for completely uniform large-scale fields in a 
triply-periodic domain. This is clearly artificial, but it provides an 
important benchmark. 

A number of surprising results have been found for weaker 

— 1/2 

fields of less than Boq. In virtually none of those cases does 

the imposed-field method recover the kinematic value of a. In- 
stead, Qimp can attain strongly suppressed values, but it can ac- 
tually also attain strongly enhanced values. This is caused by the 
unavoidable emergence of meso-scale dynamo action. In princi- 
ple, such meso-scale dynamo action could have been suppressed 
by restricting oneself to scale-separation ratios, kf/k\, of less 
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th an 2 or so. This was done, for example, in some of the runs 
of [Brandenburg &"S ubramanian (2005a). In the present case of a 
triply-periodic box, four different magnetic field configurations can 
emerge. The first three correspond to Beltrami fields, where the 
wavevector points in one of the three coordinate directions. The 
fourth possibility is also a Beltrami field, but one that varies diag- 
onally in a direction perpendicular to the direction of the imposed 
field. The latter was found to be unstable in the absence of an im- 
posed field, but they can be long-lived in the present case of an 
imposed field. 

In this paper, we have used the term meso-scale fields to re- 
fer to the Beltrami fields naturally generated by the helicity-driven 
dynamo in our system. A more general definition of meso-scale 
fields would encompass all fields that break isotropy, average to 
zero, and yet do not time-average to zero. In the absence of such 
fields, mean-field theory can be applied in a straightforward man- 
ner This is indeed the case that one is normally interested in. How- 
ever, when such meso-scale fields exist, they must be understood 
for determining turbulent transport coefficients, because those co- 
efficients apply then to the particular case of saturated meso-scale 
fields. 

The results obtained with the imposed-field method reflect 
correctly the circumstances in the nonlinear case where the a 
effect is suppressed by dynamo-generated meso-scale magnetic 
fields whose scale is smaller than that of the imposed field, but 
comparable to the scale of the domain. Especially in the case 
of closed or periodic domains the result ing a is catastrophically 
quenched, which is now well und erstood jpield & Blackmanl2002l : 
iBlackman & Brande nburg' 2002). This effect is particularly strong 
in the case where one considers volume averages, and thus ig- 
nores the effects of turbulent magnetic diffusion. With magnetic 
diffu sion included, both a and % have only a mild dependence on 
Rm l lBrandenburg et alj EoOSb). However, astrophysical dynamos 
are expected to operate in a regime where magnetic helic ity fluxes 
alleviat e catastrophic quenching; see Brandenburg & SubramanianI 
( l2005bl) for a review. 

Determining the nature of the dynamo mechanism is an impor- 
tant part in the analysis of a successful simulation showing large- 
scale field generation. Our present analysis shows that meaningful 
results for a can be obtained using either the imposed-field or the 
test-field methods provided the departure of the magnetic field from 
Bo is reset to zero to eliminate the effects of dynamo-generated 
meso-scale magnetic fields. Conversely, if such fields are not elim- 
inated, the results can still be meaningful, as demonstrated here, but 
they need to be interpreted correspondingly and bear little relation 
to the imposed field. On the other hand, for strong imposed mag- 
netic fields {RmBo/Bcq > 1), meso-scale magnetic fields tend 
not to grow, so the resetting procedure is then neither necessary nor 
would it make much of a difference when the test-field method is 
used. However, when the imposed-field method is used, the reset- 
ting of the actual field reduces the quenching of itrms. This affects 
the normalizations of Bo and aij with Bcq and ao, respectively, 
because both are proportional to itrmsH 

Throughout this paper we have considered relatively moderate 
values of Rm, but we computed a large number of different simu- 
lations. In the beginning of this study we started with larger values 
of Rm and found that the resulting aimp seemed inconsistent. In 



^ This explains why AtUrms kf is 70 in Fig. [6] and 50 in Table |4] under 
otherwise comparable conditions, except that here only the test-field is reset 
and not the actual fluctuating one. 



hindsight it is clear what happened: the few cases that we had in 
the beginning were all scattered around different branches. Only 
later, by performing a large number of simulations at smaller val- 
ues of Rm it became clear that there are indeed different branches. 
This highlights the importance of studying not just one or a few 
models of large 7?m, but rather a larger systematic set of interme- 
diate cases of moderate iim where it is possible to understand in 
detail what is going on. It will be important to continue exploring 
the regime of larger _Rm, and we hope that the new understanding 
that emerged from studying cases of moderate Rm proves useful 
in this connection. According to the results available so far, we 
can say that for larger values of R m the turbulent transport c oef- 
ficients are only weakly affected (see lBrandenburg et al.l2008bl . for 
-Rm < 600) for fie lds of equipartition strength, or not affected at 
all dSur et al.ll2008l for ii^ < 220) if the field is in the kinematic 
limit. 
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